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Abstract. We give independent proof of the deficit of stars in the in-plane central disc (2 < J? < 4 kpc, |fo| ^ 3°) 
with respect to the predictions of a pure exponential density distribution. 

We use three different methods: 1) the inversion of the red clump giant distribution in near-infrared colour- 
magnitude diagrams to obtain the star density along the line of sight; 2) the determination of the density dis- 
tribution of 1612 MHz sources by means of the distance determination of OH/IR sources from their kinematical 
information; 3) an analysis of near- and mid-infrared star counts and comparison with models. All the tests give 
the same result: a deficit of stars in the inner disc with respect to an exponential disc (either with constant 
scaleheight or extrapolated from the outer disc), but only in near plane regions ^ 3°). This deficit might be 
interpreted as a flare in the vertical distribution. The in-plane density is almost independent of R and not an 
exponential law of the type p oc exp{—R/h). Further away from the plane, however, the density increases towards 
the centre due to the increase of the scaleheight. Tests also show that this result cannot be due to extinction. 
This deficit affects both the young and the old populations, so it is probably a rather stable feature of the disc, 
and might be due to the existence of an in-plane bar sweeping the near-plane stars. 



An approximate expression of the disc density within 2 < J? < 8 kpc is: p{R, z) oc e \ « /e'^^,with 

h,{R) ^ 285[1 -f- 0.21 kpc"^ {R - Rq) + 0.056 kpc'^ {R - Ref] pc. 

Key words. Galaxy: general — Galaxy: stellar content — Galaxy: structure — Infrared: stars — radio: stars 
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K*" 1. Introduction to Anderson et al. (2002) based on analysis of the bulge- 
is^ disc decomposition carried out in Baggett et al. (1998). In 
^ Lopez-Corrcdoira et al. (2002, hereafter L02) obtained a ^ ^^^g^ ^^^^^^ external galaxies, due to their more 
^ model of the outer parts of the stellar disc (galactocentric favourable orientation, the central hole in the disc has 
radius R greater than 6 kpc), based on the 2MASS near ^^^^ detected in the distribution of molecular and atomic 
infrared point sources survey. Among their findings are the g^^^ example, Aguerri et al. (2001), who describe 
stellar warp coincident with the gas warp, the flare, the ^ ^^^^ jjj distribution in NGC 5850 associated with 
absence of an external cut-off out to i? = 15 kpc and the ^^^^^^^ ^^^-^^^ ^^^^ presented by Bottema 
measurement of the different parameters of the old pop- ^ Verheijen (2002) in NGC 3992. These holes have been 
ulation density distribution, generally m agreement with attributed to the presence of gas velocities perpendicular 
other independent studies. Here, we continue the analysis ^j^^ Galactic plane (Sancisi 1999). However, there are 
of the disc structure m the inner stellar disc, examining ^^^^ ^^^^^ ^^^^^ CO distribution, used to infer the 
whether or not this corresponds to an extrapolation of the jjjj ^^^^j^y^ ^^^^ ^^^^^1 ^^^^ ^ j^^l^^ ^^^^^ 
L02 disc model towards the centre and, if not, what the ^r example, of Sakamoto et al. (1999), whose sample is 
differences are. somewhat biased towards the non-hole gas distribution be- 
A major concern is whether extrapolation of the expo- cause of selection criteria. The hole does not usually con- 
nential in the inner stellar disc is appropriate or whether gigt strictly of an abrupt truncation of the disc. It can rep- 
there is some deficit of stars with respect to this distri- resent a gradual density decrease (for instance. Freeman 
bution. It is well known that many barred galaxies have type II truncation; Freeman 1970), instead of increasing 
holes (Ohta et al. 1990); up to half of them, according towards the center within a certain radius, or even a con- 

slant density or increasing density but with a slope much 
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lower than that produced by the extrapolated exponential 
law. 

In our Galaxy, this hole (from now on, we will call 
the hole a "deficit" of gas or stars since, strictly speak- 
ing, the hole is not totally empty; "deficit" here is always 
with respect to the extrapolation of the purely exponen- 
tial disc towards the inner radii) is observed in the gas 
distribution (i.e., Robinson et al. 1988). For the stellar dis- 
tribution the answer is not very clear at present. Several 
disc models have been produced using COBE/DIRBE flux 
maps. Although the general features are approximately 
coincident, it seems that there is no agreement regard- 
ing the deficit of stars. For instance, Freudenreich (1998) 
and Lepine & Leroy (2000) agree that there is a such a 
deficit, while Bissantz & Gerhard (2002, hereafter B02) 
can reproduce quite accurately the flux maps with an ex- 
ponential disc even in the inner parts. Caldwell & Ostriker 
(1981) gave a model of the disc that consists of a dif- 
ference of two exponentials with different scalelengths, 
which gives a deficit of stars in the inner part (see fig. 15 
in Caldwell & Ostriker 1981) based on kinematical mea- 
surements. Microlensing experiments seemed to favour a 
deficit of stars (Kiraga & Paczyhski 1994), but the poorly- 
known interstellar extinction at optical wavelengths gives 
some uncertainties in their models. Moreover, Kiraga et al. 
(1997) recognize that their previous paper might contain 
some errors due to an incorrect identification of popula- 
tions in the optical colour-magnitude diagrams. Analysis 
of OH/IR sources, both from the young and intermediate- 
old populations, led Baud et al. (1981) to conclude that 
there must be a deficit of stars in the two populations; 
these results are more convincing, and we think that a 
corroboration of this result with newer surveys will be 
valuable (see flSj. This topic has so far not received much 
attention in terms of star counts. Lopez-Corredoira et al. 
(2001, hereafter LOl) have claimed that a deficit of stars 
is necessary to reproduce some near-plane counts but they 
did not try to fit off-plane regions, and not too much at- 
tention was paid to the problem in that paper. Precisely 
because of this, we wish to examine the possible deficit 
of stars in the inner disc, using point source catalogues, 
which contain more information than flux-maps alone. 

A further area of study is the real nature of the deficit 
of stars, should this be genuine. A fiared stellar distribu- 
tion places the stars at a higher elevation from the plane 
than that predicted from a pure exponential law, hence 
giving the impression of a hole in the plane. Flared struc- 
tures are well recognized in the distribution of gas (see 
Burton & Hekkert 1986) and stars (L02; Alard 2000, and 
references therein) in the outer disc. The presence of that 
type of morphology in the central disks of spiral galaxies 
is not so well studied due to the inherent difficulties in 
separating the several structural contributions to the ob- 
served counts or flux. A good example of that difficulty 
can be found in Nikolaev et al. (2004), who reported a 
symmetric warp in the inner disk, the central 4 degrees, 
of the LMC based on accurate analysis of MACHO and 
2MASS counts. 



2. Data 

The data for this work were taken from several 
sources: mainly 2MASS, MSX and ATCA/VLA, 
and some other DENIS and TMGS data taken 
from a former paper (LOl). The second release 
of the 2MASS project (Skrutskie et al. 1997, 
http://www.ipac.caltech.edu/2mass/releases/docs.html) 
provides us with some data in the regions of interest 
to study the inner disc isolated from other components 
(1.5° < |6| < 6.5°, 15° < Z < 40°). We use these data to 
produce iiT-band star counts (mif ) and rriK vs. (J — K) 
colour-magnitude diagrams (CMDs) in the regions given 
in Tabled The MSX Point Source Catalogs (version 1.2; 
Egan et al. 1999) provides us with star counts in the D 
(14.65 /im) mid-infrared band in the regions |5| < 0.5°, 
|;| < 45°. OH/IR stars from the ATCA/VLA catalogue 
(Sevenster et al. 1997a, b, 2001) in the radio maser line 
at 1612.23 MHz at |5| < 2.5°, |;| < 45° were used to 
obtain star counts and the distribution of the sources as 
a function of galactocentric radius. 

Table 1. Selected regions used to extract the colour- 
magnitude diagram from 2MASS sources so as Ato apply 
the method of 



I (deg) 


b (deg) 


area (dcg'^) 


15 


2 


0.57 


15 
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1.00 


15 


-3 


1.00 


15 


-6 


1.00 


18 


-2 


0.92 


18 


-3 


0.97 


20 


-3 


1.00 


20 


-6 


1.00 



3. The inner disc from CMD analysis 

3.1. Method 

In L02 (§3), some of us developed a method to obtain the 
star density along a line of sight (/, b) from an analysis of 
near-infrared CMDs, such as tuk vs. (J — K), which are 
the colours to be used here. First, the trace produced in 
the CMD by the red clump giants (spectral type: K2 HI) is 
identifled. This is easy because they are much more abun- 
dant than the other giants (Cohen et al. 2000; Hammersley 
et al. 2000). Second, we count the number of stars within 
a 0.4 mag width trace as a function of the apparent mag- 
nitude (examples are given in L02 and in Figs d El of this 
paper). Third, because we know their absolute magnitude 
and intrinsic colour [Mk — —1.65, {J ~ K)o — 0.75; L02], 
we can obtain the extinction and density as a function 
of distance along the line of sight. The correction of the 
mean extinction is carried out since we can measure the 
reddening of the K2 III stars. The details are explained in 
L02(§3), so they will not be repeated here. 
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Recent results in open clusters also give similar results 
for the intrinsic luminosity of the red clump (Grocholski 
& Sarajedini 2002): magnitude —1.62 with a standard de- 
viation of 0.21 and colour (J — K)q around 0.7 for the 
age of the old disc population with a small dispersion due 
to metallicity or age gradients. Pietrzynski et al. (2003) 
also show that the i^-band magnitude of red clump is 
a good distance indicator, and that the colour has only 
a very slight dependence on the metallicity. Theoretical 
predictions from isochrones (Salaris & Girardi 2002) agree 
with these values. Therefore, counts can give us informa- 
tion about the extinction along the line of sight as well 
as the density of K2III stars: pK2Ui{i",l,b), which leads 
to Pk2IIi(^j^)j assuming axisymmetry in the disc (r = 
distance from the Sun, R = distance from the centre of 
the Galaxy). 

Concerning the completeness in crowded fields, in the 
worst cases {I = 15°, b — 2°) the survey is complete to 
magnitudes ttik ~ 13.2 and mj w 14.4. The reddening 
due to the extinction is less than A(J — K) — 0.75. Since 
we use the red clumps with (J — K)o between 0.55 and 
0.95 (the width of the strip is 0.4 mag), we will adopt a 
limit of rriK < 12.7 (instead of 13.0 used in L02) to be 
sure of the completeness of the sources. 

Poisson errors will be important for very bright stars 
(i.e. close stars, since the absolute magnitude is constant). 
If we look towards the centre of the Galaxy, the Poisson 
error will be important for distances from the Sun r < 3 
kpc (galactocentric distance i? <; 5 kpc), but it will be 
reasonably small for R — 2-5 kpc, where the inner disc 
can be examined. The method works reasonably well for r 
greater than ~ 3 kpc, and this was tested with real data by 
L02, who demonstrated that the parameters of the outer 
disc obtained by this method are consistent with those 
obtained by other methods. The global systematic errors 
of this method are less than 10% for the regions with low 
extinction in a typical anticentre direction (L02). In the 
present paper, the areas are somewhat more extinguished 
(up to 0.45 mag in K instead of 0.2 mag in L02, total 
extinction out to 7-8 kpc along the line of sight) and the 
regions are towards the centre of the Galaxy instead of the 
outer disc regions used in L02, which might give a some- 
what higher contamination of M-giants. This, however, 
will give systematic errors of the same order, as will be 
shown below, perhaps only slightly larger, and the method 
is therefore appropriate for measuring the density of the 
inner disc stars. 



3.1.1. A simulation to test the method 

To show how well the method works, we will carry out 
a Monte Carlo simulation. We will generate a synthetic 
colour-magnitude diagram with a population model and 
a predescribed density model. From this CMD, we will 
obtain the density along the line of sight using the method 
of extracting the clump giants, and this will be compared 
with the original density. 
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Fig. 1. Top: Colour-magnitude diagram for the line of 
sight I = 15°, & = 2° used in the present section with 
2MASS data. The solid line shows the fitted trace that we 
assign to the red clump giant population and the limits for 
the red clump giants within a width of 0.4 mag. Dashed 
lines show the limits for a width of 0.3 mag. Dot-dashed 
line show the limits for a width of 0.6 mag. Bottom: his- 
tograms of counts (per unit colour) for three fixed appar- 
ent magnitudes and the corresponding best Gaussian fits: 
rriK— 11-95 (dots/solid hnc), mx=12.25 (crosses/dashed 
line), rriK — 12.55 circles/dot-dashed line. Am = 0.3. 

The information about the absolute magnitudes and 
relative abundances of all possible stellar spectral types 
was taken from the updated SKY model (Wainscoat et al. 
1992; M. Cohen, private communication). The red clump 
with the range of intrinsic colours of 0.2 mag and range 
of absolute magnitudes of 0.3 mag (L02) is represented 
by the K0-5III population with a peak in K2-3III; this 
population represents the clump with mean absolute mag- 
nitude — 1.65 and color (J — K)q — 0.75, with Gaussian 
dispersions with the referred sigmas. Poisson noise in the 
density is also introduced in the simulation. 

Since the region at Z = 15°, 6 = 2° has a higher 
contamination of M-giants (because the ratio of dis- 
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Fig. 2. Top: Colour-magnitude diagram for the hue of 
sight I = 20°, 6 = — 6° used in the present section with 
2MASS data. The solid line shows the fitted trace that we 
assign to the red clump giant population and the limits for 
the red clump giants within a width of 0.4 mag. Dashed 
lines show the limits for a width of 0.3 mag. Dot-dashed 
line show the limits for a width of 0.6 mag. Bottom: his- 
tograms of counts (per unit colour) for three fixed appar- 
ent magnitudes and the corresponding best Gaussian fits: 
rriK = 11.95 (dots/solid line), niK — 12.25 (crosses/dashed 
line), mx=12.55 circles/dot-dashed line. Am = 0.3. 

tant/nearby stars is larger) and the highest extinction, 
we will carry out our experiments in this line of sight. 
This is the one with the largest errors, so the other lines 
of sight will have lower errors than the present one. We 
carry out two experiments: 1) with a pure exponential 
disc: piR,z) cx e-^/^e"!^!/''- (i?© = 7.9 kpc, scale- 
length H = 0.25i?o, scaleheight = O.O36i?0; L02); 2) 
for a model with the same density for R > A kpc, but 
p{R,z) oc e~l^l/''^ (i.e. independent of galactocentric dis- 
tance) for R < A kpc. The results are shown in Figs|21and 

El 

The colour-magnitude diagrams in Figs |3| and 0] ap- 
proximately reproduce the one observed in Fig. ^ They 




Fig. 3. Top: Synthetic colour-magnitude diagram derived 
from the SKY model for the populations and relative 
abundance and an exponential disc (larger density for low 
values of R). Centre: histograms of counts for fixed mn 
(Amx — 0.3). Bottom: density along the line of sight in 
the field I = 15°, & = 2°: comparison between the "real" 
assumed density and the corresponding densities obtained 
from the CMD analysis method with width = 0.4 mag of 
the synthetic colour-magnitude diagram. 



are not perfectly equal, but it is not our purpose to discuss 
the validity of the SKY model here. Our purpose is to test 
the density extracted through the clump giant extraction 
method. This is tested in the right-hand part of Figs 01 
and El (for r > 3 kpc; for lower distances, the Poisson er- 
ror and the M-giant contamination dominate, but we are 
not interested in this range). In both cases, the density 
is approximately that of the model. There are excesses of 
20-30% between r — 4-6 kpc due the contamination of 
the M-giants. The density distributions in both cases can 
be clearly distinguished and the method is able to detect 
the differences of the density in Figs O and El This is the 
worst case among the selected lines of sight. Statistically 
(when we use the data of all the lines of sight), the errors 
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"SKY" + L02-constant dens, for R<4 kpc 




Fig. 4. Top: Synthetic colour-magnitude diagram derived 
from the SKY model for the populations and relative 
abundance, and an exponential disc for i? > 4 kpc and 
constant density (independent of R) for i? < 4 kpc. 
Centre: histograms of counts for fixed mi^ {Anij^ — 0.3). 
Bottom: density along the line of sight in the field I = 15°, 
6 = 2°: comparison between the "real" assumed density 
and the corresponding densities obtained from the CMD 
analysis method with width — 0.4 mag of the synthetic 
colour-magnitude diagram. 



are lower: roughly 10-15% for r 
beyond. 



4-6 kpc and even less 



3.1.2. Uncertainties due to the patchiness of extinction 

The extinction is corrected on average because we know 
the reddening. In extracting the red clumps from the 
colour-magnitude diagram, we ask whether the patchi- 
ness of the extinction, i.e. the fluctuations in relative 
to the mean extinction, could have some effect (charac- 
terized by dispersion (tak)- This is an additional effect 
to consider which was not discussed in H3.1.1I It would 
produce an extra dispersion in the apparent magnitudes 
(which, if less than ~ 0.3 mag, produces negligible er- 



rors in an exponential distribution; see L02, §3.3.1) and 
an extra dispersion in the reddening, which can result in 
the loss of some stars because they move to a position in 
the CM diagram that is outside the selected trace of red 
clump giants with width 0.4 mag (see L02, §3). Since the 
intrinsic dispersion of colours is around 0.2 mag, the extra 
dispersion with respect to the central position of the trace 
would be (for A(J - K) = 1.52^^^; Rieke & Lebovsky 
1985): 



f(J--fs:), centre trace — ^-52(7 Ak 



( drriK \ 
\d(J-K)J 



(1) 



where the first term is the absolute value of the reddening 
and the second term is the (J — K) variation of the centre 
of the trace corresponding to this extra extinction, j^^^"^^ 
is the equation of the centre of the trace. From the relation 
between apparent magnitude and the extinction, we get 



dm 



K 



d{J -K) 



drriT. 



d(.l-K) 
dr 



Slogio e 



dAK 
dr 



1.52 



dAi. 



= 0.658+- 



1.43 

dAK(r) 
dr 



(2) 



Consequently, for instance for r = 5 kpc at / = 15°, h = 

2° (corresponding to i? = 3.4 kpc): Cr(j_ if), centre trace = 

l.liaAK- For larger values of r (lower R) the effect is 
somewhat lower because the second term of eq. be- 
comes more important: for instance, cr(j_i<-), centre trace — 
0.92(7^^ (for r = 7 kpc). This extra broadening of the dis- 
tribution will lead to a loss of stars for a fixed width of the 
trace. Assuming that this extra broadening and the intrin- 
sic broadening of the distribution are Gaussian, the ratio 
of lost stars would be that plotted in Fig. Ofor I ~ 15°, 
b = 2°. 

A Monte Carlo simulation was also carried for the cases 
of Figs 121 and 21 but adding the effect of patchiness of ex- 
tinction. Fig. (HI shows the result. The effect of the patch- 
iness of extinction might be important, but only in cases 
with uak = 0.6Ai<- or larger, while the expected cases are 
under Q.2Ak, as will be shown now. 

To determine whether the patchiness of extinction is 
significant it is necessary to evaluate the dispersion of this 
extinction. An examination of the Schultheis et al. (1999) 
and LOl maps of extinction (with a resolution of 2 arcmin- 
utes) for |Z| < 20°, |6| < 2° in Ay gives a spread (due to 
both patchiness and gradient) of around 20% with respect 
to the mean extinction, almost independently of which re- 
gion in the centre of the Galaxy is considered. 

Lada et al. (1999) have shown that the variations 
in extinction through the dark cloud IC 5146 for each 
background star due to small-scale fluctuations (once the 
large-scale gradients are removed, i.e. taking small re- 
gions) are less than 18% of the mean extinction. If this is 
correct for the Milky Way extinction in general, it would 
indicate that the patchiness is limited within the value of 
c^Ak — 0.2QAk, or even less if we integrate along a line of 
sight that passes through several clouds, as would be the 
case for near-plane regions. 
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Fig. 5. Fraction of lost stars due to the spread of the ex- 
tinction in the field / = 15°, b = 2° for r = 5 kpc. Two 
values for the width of the K2III strip are used. In the 
subsequent analysis, a width of 0.4 mag is used. 



Since our region has 0.5 mag of mean extinction, this 
would give ~ 0-10 mag, i.e. a loss of 10% of the stars 
in our case (7% for width 0.6 mag). This would be nearly 
the same for any r between 5 and 6.5 kpc (i? between 2.3 
and 3.5 kpc). This is the worst of the cases, because the 
rest of the lines of sight are less extinguished and are con- 
sequently less patchy. A model with a substantial deficit 
of stars (such as the one we will obtain below) gives an 
absence of 60-70% of stars with respect to the purely expo- 
nential L02 model at this latitude. A value of around 
0.40 mag, four times our estimate, would be necessary to 
explain this, which makes this method able to detect a 
deficit of stars of this order. 

The effect of the gradient of extinction may be shown 
to be insignificant overall. For the most extinguished re- 
gion {I — 15°, b = 2°), the effect of the gradient of ex- 
tinction due to variation of b between 1.5° and 2.5° is 
^ = Ab^^MlL^ and -^^'"/^^ « 0.49 (Wainscoat et 
al. 1992); the gradient in extinction in the latitude di- 
rection produces a dispersion of ~ 0.07 mag, which 



is certainly not Gaussian but for a rough estimate can 
be considered so. From Fig. [S] we will see that this can 
produce a 5-7% loss of stars, which may be considered 
unimportant. 

3.2. Application of the mettiod 

We selected eight regions (see Table in |^| = 15-20°. 
In this zone, the disc is far enough from the bulge that 
the contribution of bulge stars in the extracted red clump 
giant strip is small, and it is close enough to the Galactic 
Centre to have lines of sight reaching R — 2-4 kpc, where 
the inner disc will be studied. 

Different latitudes were selected so as to have a rep- 
resentation of different heights, except in plane regions, 
where the extinction might be excessive and/or very irreg- 
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Fig. 6. Density along the line of sight in the field I = 15°, 
b — 2° [a) exponential disc as in Fig. O b) with con- 
stant density for R < 4 kpc as in Fig. 0]: comparison 
between the inputs of the density and the density distribu- 
tion obtained from the CMD analysis method with width 
= 0.4 mag of the colour-magnitude diagram generated 
by a Monte Carlo simulation. Models with four different 
patchiness are plotted: cr^^. ~ 0, 0.2^^, 0.4A/f , 0.6v4^. 



ular in spatial distribution, and the hypothetical in-plane 
bar would presumably contaminate the counts (LOl). 

The density pk2Ui{Rt z) along the eight lines of sight 
was obtained (see Figs [TJ . An extrapolation of the ex- 
ponential L02 model does not fit these data: the model 
gives only '^30% of the density at 6 = —6°, which means 
that the extrapolation of the flare (the decrease in with 
decreasing R) found in L02 for the external disc is not 
appropriate for the inner disc. If we extrapolate the expo- 
nential L02 disc but without a flare, the result is better 
for b — —6°, but there is an excess of ^200% in the den- 
sity at |6| = 2°. It is clear that an exponential disc does 
not work either with or without a flare of increasing h^ 
for increasing R, according to the extrapolation of L02, or 
with any intermediate option. 
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Fig. 7. Density along the different lines of sight. An ex- 
trapolation of the exponential L02 model does not fit 
these data: the model gives only ^30% of the density at 
6 = —6°, which means that the extrapolation of the flare 
(the decrease in hz with decreasing R) found in L02 for 
the external disc is not appropriate for the inner disc. To 
fit these data, a density nearly independent of R must be 
introduced (thicker solid line in the plot), which will be 
obtained in this paper [eqs. Q and ©]. 



In Fig.|Sl the densities are plotted as a function of z for 
different R in R < 4 kpc (note that, for the calculation of 
z(r, I, b), we take into account that the Sun is 15 pc above 
the plane, but not the effects of the warp described in L02 
since they are negligible for the inner disc; the warp is 
important for regions near |/| = 90° and i? > 10 kpc). In 



these plots, we observe how the scaleheight changes with 
R; it increases towards the centre: hz{2.5 kpc) = 424 ±46 
pc, hzi3 kpc) = 352 ± 23 pc, /i^(3.5 kpc) = 338 ± 48 pc, 
/iz(4 kpc) = 323 ± 26 pc. A weighted fit of a linear law of 
hz gives: 



317 ± 17 - [i?(kpc) - 4]48 ± 20 pc. 



(3) 



for 2.25 kpc < R < 4.25 kpc. Therefore, we have a higher 
scaleheight in the inner disc with the opposite trend to 
that of the outer disc, which, according to L02, ranges 
from 230 and 285 pc between R — 6 kpc and Rq: here, 
hz slightly decreases with increasing R. Eq. Q is valid 
for R < 4.25 kpc while that given by L02 is valid for 
R > 6 kpc. In the range between 4 and 6 kpc the scale- 
height dependence on radius is somewhat uncertain but is 
presumably the region where hz is at a minimum. 

One might suspect the contamination of the bulge 
to be responsible for this increase in the scaleheight. 
However, the bulge contribution is negligible. Using the 
bulge model of axial ratios 1:0.54:0.33 of Lopez-Corredoira 
et al. (2000), we find that the largest bulge/disc ratio is 
for I = 15° , b ~ —6° , R = 2245 pc which gives a contri- 
bution pbuigc = 5 X 10^"^ pc~'^. Since the K2 giants are 
sa 10~^ times the total number of stars (L02), this gives 
/3K2iiLbuigc ~ 5 X 10^^ pc^^ = 5% of thc total K2III 
density, which is less than the other errors. Higher ellip- 
ticity models of the triaxial bulge (such as that of B02) 
would give even less bulge contamination. This is the most 
favourable line of sight for the bulge; for the other (/, b, R) 
the bulge contamination is even lower. 

Once we know the scaleheight as a function of the ra- 
dius, hz{R), we can calculate 



PK2m[R,z = 0) = -rji]- , 



(4) 



gh^(H) 



where hz (R) is taken from eq. for i? < 5.9 kpc and from 
L02 for R > 5.9 kpc, in order to keep the continuity of hz- 
This is represented in Fig. |5| For i? > 5 kpc the dispersion 
is larger because the points stem from the counts of very 
bright red clump sources (which are few in number), so 
they have significant Poisson noise and also an important 
contamination of M-giants. It is remarkable how low the 
dispersion of points is for i? < 4 kpc and how clearly 
the density does not increase exponentially towards the 
centre. If eq. representative of the scaleheight in the 
inner disc, was greatly in error, the dispersion for i? < 4 
kpc would be much higher, since we have used data with 
|5| between 2° and 6°. 

If we fit a linear relationship in Fig. |51 for R < 4.25 
kpc, we obtain: 



PK2iii(i?,0) = (4.57±0.33)[1 
- 0.150 ± 0.033(i?(kpc) - 4)]pK2iii,0; 



(5) 



for 2.25 kpc < R < 4.25 kpc, where pk2IILo = 1-31 x 
10"^ pc~'^ (L02). This shallow rise is incompatible with 
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Fig. 8. Density of red clump giants, pk2IIi(-R, z), at R —2.5, 3.0, 3.5 and 4.0 kpc. The intervals of R are in steps of 
0.5 kpc, and the data in each diagram belong to i? — 0.25 kpc < galactocentric radius< R + 0.25 kpc. Circles: line of 
sight I = 15°, b = 2°; squares: / = 15°, b = 3°; diamonds: I = 15°, b = —3°; triangles up: I = 15°, b = —6°; triangles 
left: I = 18°, b = -2°; triangle down: I = 18°, b = -3°; triangle right: I = 20°, b = -3°; plus: I = 20°, b = -6°. Best 
exponential fit vs. z is also plotted. 



an extrapolation of the outer exponential disc towards the 
centre. Although the data to 2 kpc show a slight increase of 
the density towards the centre, a constant density is also 
consistent within the systematic errors {Ap/p ^ 20%). 
If we neglect the small gradient with the galactocentric 
radius in the plane (since the systematic errors are larger 
than these small variations and we cannot be sure that 
the gradients are real), then 

PK2IIi(-R, Z = 0) « 5.1pK2IILO, (6) 

for 2.25 kpc < R < 4.25 kpc. We can also fit a more 
complex law, like a modified exponential (already used 
for instance by Lepine & Leroy 2000): 

p(i?,z = 0) =pQeV « +«0''e"i«+— (7) 



where pq — 1.31 x 10^^ pc~^ and H — 1.97 kpc are 
already known from L02, and the parameter H2 is fitted 
with our data (see Fig. O with a best value of: 



H2 = 3740 ± 130 pc. (8) 

The fit of eqs (0) and (0) , or eq. © , together with eq. 
(PJ, is also applicable away from the plane, as can be ob- 
served in Fig. ^] Observe that for lower z, the observed 
density is lower than the extrapolation of the exponen- 
tial law, while for higher z the density increases towards 
the centre more rapidly than the extrapolation of the L02 
model. That is, the deficit of stars with respect to the L02 
disc affects the regions near the plane. Also, Fig. [Tlshows 
that the present model fits the density obtained directly 
from the counts (as a function of the distance from the 
Sun for the different lines of sight). 
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Fig. 9. Density of red clump giants in the plane [derived 
indirectly from the counts, using the scaleheight hz{R) 
from eq. Q for R < 5.9 kpc and from L02 for R > 5.9 
kpc]. Note that the data for R > 5 kpc (r < 3 kpc) have 
large errors due to Poisson noise and M-giant contamina- 
tion. Exponential disc model taken from an extrapolation 
of L02 outer disc model. Circles: line of sight I = 15°, 
b = 2°; squares: I = 15°, b = 3°; diamonds: I = 15°, 
b = —3°; triangles up: I — 15°, b — —6°; triangles left: 
I = 18°, 6 = -2°; triangle down: I = 18°, b = -3°; trian- 
gle right: I = 20°, b = -3°; plus: / = 20°, h = -6°. 
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3.2.1. Flare in the inner disc 
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Another way to interpret the observed effect is the exis- 
tence of a flare in the inner disc, which redistributes the 
stars of the plane in off-plane regions. In L02 we consid- 
ered an exponential distribution for the density where the 
flare affects both the scaleheight (which increases with the 
galactoccntric distance) and the flare scalelength (which 
changes, causing the scalelength of the disc space den- 
sity, H, to remain constant; see L02 for details). However, 
the observed effect in the in-plane density at the central 
Galaxy can also be reproduced by an inner flared disc. If 
we consider a disc with a fixed exponential scalelength of 
the disc surface density, and we flare this disc vertically by 
the flaring law defined in eq. Q, conserving the mass, i.e., 
the stellar surface density, we find a radial law for the disc 
density consistent with the observed counts. Considering a 
fixed scalelength of hji — 2.4 kpc, for example, the agree- 
ment is quite good between the observed and predicted 
densities in the plane in the range 2-8 kpc. As a compar- 
ison, in L02 the scalelength of the disc surface density is 
Hr — 3.2 kpc in the outermost disc {R >10 kpc), but the 
flare scalelength changes as the intrinsic scalelength of the 
disc varies. 

In Fig. 121 we represent the predicted densities for four 
different models: i) the extrapolation of the exponential 
L02 model, ii) a linear relationship for R < 4.25 kpc as 
given by eq. (0) , iii) the modified exponential given by eq. 



Fig. 10. Density of red clump giants at fixed z = 150 pc 
and 400 pc [derived indirectly from the counts using the 
scaleheight hz{R) from eq. Q for i? < 5.9 kpc and from 
L02 for R > 5.9 kpc]. Symbols are the same than in Fig. 

El 

and iv) a flared disc with a fixed scalelength of the 
surface density (h^ = constant). In each case, except for 
the extrapolation of the L02 model, there is a good match 
between the predicted and observed density. This agree- 
ment is also observed away from the plane, where models 
ii), iii), and iv) predict roughly the same star density (Fig. 

iini. 

However, there are basic differences between a vertical 
flare and a "hole" in the central disc, the latter being con- 
sidered as a deficit in the stellar density. While the former 
retains the vertical column density, in star numbers, a hole 
in the sense defined here will decrease that column den- 
sity, and will appear to an external observer looking at the 
Galaxy face on as a real deficit in the stellar distribution 
of the central regions as compared to the surroundings. 

We plan to investigate further the real nature of this 
feature of the stellar distribution in a subsequent paper. 
In the rest of this paper, we will refer to this observed 
feature as a "hole" , without prejudging the true origin of 
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this decrease in the stehar distribution with respect to the 
exponential prediction. 

3.2.2. Calculations with different band widths in the 
subtraction of red clump giants 

The variation in the width of the band to extract the red 
clumps could alter the results slightly due to the intro- 
duction of certain systematic errors. In this subsection, 
we will examine how large these effects are. 

We have already seen f tl3.I.I(l that the contamination 
of the giants cannot emulate such a strong deficit of stars. 
Furthermore, the contamination would increase the counts 
with respect the number expected in the L02 model (and 
the normalization of this model was carried out in re- 
gions where the contamination is very low) but never de- 
crease them. A more important source of error could be 
the higher loss of stars for fainter magnitudes due to the 
spread away from the 0.4 mag width strip, but we have al- 
ready verified fi l3.1.2|l that this effect is unimportant (less 
than 10% according to calculations). 

For a further test of the results, we have checked that 
changing the aperture of the trace (which also changes the 
contamination of other sources: dwarfs and M giants; see 
L02, §3) does not significantly change the density obtained 
for the inner disc. We have derived expressions ||2J) and (|SJ) 
for width 0.6 instead of 0.4 used above, with the result: 

/i^ = 282 ± 9 - [i?(kpc) - 4]48 ± 15 pc; (9) 

PK2iii(i?,0) = (5.55±0.34)[1 

-0.024 ± 0.020(i?(kpc) - 4)]pK2iii,0; 

for 2.25 kpc < R < 4.25 kpc, (10) 

which is compatible with the previous result. If we re- 
duce the width rather than increasing it, for instance, to 
a 0.3 mag width instead of 0.4 mag, the results are: 

= 325 ± 16 - [i?(kpc) - 4]49 ± 21 pc; (11) 

PK2iii(i?,0) = (5.15±0.40)[1 

-0.165 ± 0.037(ii(kpc) - 4)]pK2 111,0; 

for 2.25 kpc < R < 4.25 kpc. (12) 

There are some small differences that can stem from 
the different degree of contamination from stellar types 
other than K2III (a width of 0.6 is more contaminated), 
and the different degree of patchiness of extinction (this 
affects the width 0.3 more; see ii3. 1.2(1 . These small differ- 
ences do not change any of the conclusions. 

On the other hand, the possible errors in the char- 
acterization of the K2III stars would move the position 
of the beginning of this non-exponential distribution, but 
never reproduce the shape of pk2IIi(-R, z = 0) decreasing 
towards the centre. 



3.2.3. Summary of the results of this section 

A rough expression for the central inner disc density is: 

p{R,z) w 5.1p0e(™9-4SBrkpc) PC) (13) 

This is applicable to 2.25 kpc < R < 4.25 kpc, and the re- 
gion of transition with the outer L02 disc would be abrupt. 
A more realistic model is likely to have a smooth transi- 
tion between 4 and 6 kpc, perhaps something like eq. ((T^ . 
This expression, eq. (|13|l . is equivalent to: 

-(R~R0) -w 285 pc 

p{R,z)K.pQe 2.4 kpc e(509-48fl(kpc) PC) ^ (14) 

hz{R) 

a flared disc with a fixed scalelength of 2.4 kpc and a 
scaleheight that varies according to eq. 

Due either to a central hole in the stellar distribu- 
tion or to a flared disc that redistributes the stars in 
higher heights above the plane, there is a quite signifi- 
cant deflcit of K2 III stars in Galactic plane in the inner 
disc. Therefore, although we must treat these results with 
some care, we think that the present conclusions of this 
section are to be considered as a tentative detection of a 
flat (rather than exponential) density distribution in the 
inner disc of our Galaxy for near plane regions. A remark- 
able aspect of this deflcit of stars is that it affects mainly 
low z, i.e. low-latitude regions, because the larger scale- 
height in the central regions compensates for the deficit of 
star for higher latitudes (see Fig. llOfl . This agrees with the 
fact that the deficit was not detected in Baade's windows 
(at h — —3.9°) (Kiraga et al. 1997). High z regions in the 
model do not have a large deficit of K2 III stars but even 
an excess with respect to the extrapolation of L02 outer 
disc model (see Fig. ^] for z — 400 pc). Presumably, the 
deficit of stars affects only the near-plane regions, where 
the disc material can be swept up due to the motion of 
the in-plane bar (Athanassoula 1992, LOl). 

4. Star counts in the near-infrared 

A fiat density distribution in the inner disc would also 
produce a deficit of star counts relative to an exponen- 
tial disc. Counts up to magnitude 9 in the K band are 
representative of the old stellar population and are ap- 
propriate for studying the central regions of the Galaxy. 
Much deeper counts would make the local disc predomi- 
nant instead (Garzon et al. 1993). 

In Figs E] and El we see how the model with a fiat 
density distribution fits the counts in off-plane regions 
quite well, while the extrapolation of the L02 external 
disc model or the B02 disc model (the luminosity func- 
tion and the extinction were taken from the SKY model; 
updated version of Wainscoat et al. 1992) is significantly 
worse, especially for near-plane regions. All models were 
normalized to give the best fit to the counts presented in 
Fig. ^] The extinction is very small, so it does not matter 
whether the applied corrections are totally accurate. No 
attempt is made to fit the plane {z = 0) regions, where 
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Fig. 11. 2MASS star counts with < 9.0 for 15° < I < 
40°. Comparison with B02 disc model, and the extrapo- 
lation of the L02 model of outer disc towards the centre 
with exponential and flat density distribution (this paper) 
are shown. All models were normalized to give the best fit 
to the counts. 

the in-plane bar (LOl) and the higher extinction make the 
conclusions about the disc less clear. 

Since the L02 model is based in star counts, it rep- 
resents better than B02 (based on flux maps) the varia- 
tion of counts with latitude. L02 model also includes the 
variation of hz{R) (a flare), making hz lower towards the 
centre, improving the fit. However, the L02 extrapolation 
is insufficient to fit the counts, especially in 6 = —1.75° 
and also quite remarkably in 6 = —4°; for higher latitudes 
there are not so many differences. From here it seems clear 
that both the flare and the inner truncation are necessary 
to fit the counts. 

Obviously, in Fig. ^] we see a high excess of counts 
in 1^1 < 10° which is not reproduced by the disc model 
because it is caused by the bulge. We see also some excess 
in positive longitudes with respect to negative longitudes, 
probably due again to the in-plane bar (LOl). The most 
remarkable fact about the disc in Fig. E|is, as was said 
in LOl, that the counts are constant between |/| = 15° 
and \l\ — 30° (it is clearer at negative longitudes), and 
this cannot be fitted with a purely exponential disc in the 
inner stellar disc. COBE/DIRBE data at 2.2 /im also give 
a nearly flat light distribution, which at least does not 
decrease towards the centre, for 15° < |/| < 30°, 6 = 0° 
(Hammersley et al. 1994, Fig. 1); although in 6 = 0° the 
imprint of other features such as the in-plane bar or spiral 
arms also appears. 

4.1. Discussion of extinction effects 

For the near infrared counts, particularly for |6| = 1.75°, 
where the deficit of counts in the central part is more evi- 
dent (see Fig. ll2() . the extinction is not very high, making 
it rather unlikely that this is a problem of the calcula- 
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Fig. 12. DENIS/TMGS star counts with < 9.0 for 
-2° < 6 < -1.5° (LOl). Comparison with the B02 disc 
model, and the extrapolation of the L02 model of outer 
disc towards the centre with an exponential and flat den- 
sity distribution (this paper) are shown, with the same 
normalization as in Fig. ^2 Note that the models repre- 
sent only the disc (without the bulge, bar and ring) while 
the counts do include everything in the line of sight. 



tion of the extinction. For instance, at |/| = 15° there is 
a deficit of 30-45% in the counts with respect to the L02 
model without a flat inner density (and even more for the 
B02 model), which would need an error of 0.3-0.4 mag in 
extinction in order to be explained in terms of dust. With 
the model that we have used, the total extinction in the 
Galactic inner disc at I = 15°, 6 = 1.75° is 0.48 mag (up 
to a distance from the sun of 7.6 kpc). This is consistent 
with the estimate of the reddening for the red clump stars 
in the colour-magnitude diagrams (see ^ at I = 15°, 
6 = 2°, which gives wO.45 mag of total extinction up to 
that distance (for 6 = 1.75° it would be ^^0.06 mag larger). 
Another independent estimate based on Schultheis et al. 
(1999) gives Ay « 5 — 6 mag for regions in 6 = —1.75°, 
7.5° < / < 17.5° (see Fig. 8 in LOl), i.e. 0.5-0.6 mag in K, 
smaller than the ~0.8-0.9 mag or even more, necessary to 
produce the observed deflcit of stars. A further difficult 
point, in case we suspected that the extinction is respon- 
sible for this distribution of counts, would be to explain 
the constant star counts in 15° < \l\ < 30°. It is difficult 
to explain an excess of extinction precisely in a region 
where we know that there is a hole of gas/dust along the 
line of sight. The possible extra extinction produced by a 
ring of gas of radius ~3-4 kpc would give less extinction 
for lower |/|, because the crossing of the line of sight in 
the ring is more perpendicular to the ring than tangen- 
tial; and, moreover, the ring is likely to have a significant 
contribution only for |6| < 1° (LOl). Therefore, we feel 
that many coincidences that would be difficult to explain 
would need to take place in order to explain the observed 
facts in terms of only a problem of extinction calculation. 
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Fig. 13. DENIS star counts with TTiK.corroctod < 9.0 and 
(J — K) > 0.5 based on LOl-Fig. 10 (see method of ex- 
tinction correction in LOl, using {J — K) colors). Al = 2°, 
Ab = 0.5°. 

Examination of the counts corrected for extinction is 
a further test that can be carried out. In Fig. 10 of LOl, 
we can see the counts for the negative latitudes after they 
were approximately corrected for extinction (see LOl for 
details). The colour (J — K) > 0.5 was chosen to remove 
part of the local disc sources and improve the contrast 
of distant sources. In Fig. E|we reproduce the counts of 
LOl (fig. 10) as a function of longitude for four different 
latitudes. In this figure, we observe first of all that the 
stellar 3 kpc arm (that is, the ring, at I « —22°) affects 
only those counts for \b\ < 1°. Second, and more impor- 
tantly, we can see that the counts are nearly longitude- 
independent for the range —32° < I < —16° at 6 = —1.5°. 
This is also true for latitudes closer to the plane, but the 
3 kpc arm at I « —22° is added. Therefore, we see again 
that the shape of the counts is incompatible with an expo- 
nential disc without truncation, which would give increas- 
ing counts towards the centre; in this case, we are talk- 
ing about counts corrected for extinction. Again, many 
unlikely coincidences would need to occur to explain the 
observed facts without the disc truncation. 

5. Density distribution of OH/IR stars at 1612 
MHz through kinematical distance 
determination 

Stars with high mass-loss rates produce OH masers whose 
strongest line is at 1612.23 MHz (Elitzur et al. 1976). 
This allows us to observe these stars at this radio fre- 
quency with the advantage of negligible extinction. A 
complete survey of OH/IR stars in the region |/| < 45°, 
|6| < 3.25°, ATCA/VLA OH 1612 MHz, was done by 
Sevenster et al. (1997a, 1997b, 2001). These sources trace 
the intermediate-old population when the sources with 
wind outflow velocity (ug) less than 14 km s^^ are se- 
lected, as will be done in this paper; they have an age 
of » 1.5 — 8 Gyr (Sevenster 1999). The completeness is 



very close to 100% for sources with flux larger than 0.375 
Jy (Sevenster et al. 1997a); therefore, we use only sources 
above this limit in this paper. 

Star counts cannot be used as a test to discriminate 
among the models since the number of sources is very 
small and the Poisson errors dominate. Another direct way 
to check whether there is a deficit of sources at i? < 4 kpc 
is to obtain the distances of the sources. The survey has 
information on the radial velocities of the sources so, in 
principle, it is possible to carry out such a test. The only 
difficulty is that the number of sources is not very high, 
although high enough to discriminate among the differ- 
ent models. The test is performed as follows: first we cal- 
culate the galactocentric distance of each source in cer- 
tain regions of the sky (which contain inner disc sources 
but avoiding regions where other components contribute a 
large number of sources); second, we count the number of 
sources in the inner disc region in successive radial rings; 
third, we compare the results with the predictions of dif- 
ferent density models using a luminosity function that we 
also obtain directly from the data. 

5.1. Galactocentric distance determination 

The first step, the galactocentric distance (i?) determi- 
nation, is carried out through the well known relation 
between radial velocity, Vr, and R for in-plane regions 
(Burton 1988, eq. 7.4), assuming circular orbits: 



R 



Rr. 



sinL 



(15) 



where V{R) is the rotation velocity of the disc at galacto- 
centric distance R. If we take V{R) = V{RQ){R/R0f-^ 
km s""'^ (Binney et al. 1991), this leads to 



R = 



Rr. 



1 



V{RQ)sinl 



(16) 



This equation is applicable because the disc OH/IR stars 
in the plane follow the galactic rotation (Baud et al. 1981). 
We adopt V{Rq) = 200 km s'^ (Honma & Sofue 1996) 
(a value of 220 km s~^ would not change the conclusions 
of this section, but would slightly reduce the values in 
Fig. Uni/ Table 13) and Rq = 7.9 kpc (Lopez-Corredoira et 
al. 2000). The random motions in the old population are 
around 30-40 km s'^ (Baud et al. 1981). cr„^ « 35 km s'^ 
leads to an error of 



0.020 

sin I 



i?(kpc)^-y kpc. 



(17) 



and this leads to important errors only for large values of R 
and low I. Moreover, the errors for the individual sources 
are mostly reduced when the statistical count is carried 
out, except perhaps for the very few sources close to the 
Sun (high R) . In any case, we will include the calculation 
of this error in this paper to show that even with this error 
(and the Poisson errors) models with a purely exponential 
disc are excluded. 
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5.2. Luminosity function 

The previous relation allows us to determine R, but not 
the distance from the Sun (r) because of the well-known 
ambiguity of there being two possible values for a given R 
inside the solar circle (Burton 1988). Therefore, we cannot 
use this relation to obtain the luminosity function directly. 

To estimate the luminosity function, we take sepa- 
rately the sources in the centre of the Galaxy (|Z| < 3°, 
l&l < 3°) and assume that all of them are at the same dis- 
tance (i?o = 7.9 kpc). The contribution of foreground disc 
sources is small. To reduce this disc contamination fur- 
ther, we take the sources with radial velocities \vr\ > 30 
km s~^, since the disc sources in this direction are ex- 
pected to have near-zero radial velocities. Fig. El shows 
the cumulative luminosity distribution of this subsample 
(luminosity in units Jy kpc^, i.e. the flux of the source in 
Jy considering the source located at a distance of 1 kpc) 
which can be fitted by: 



ATCA/VLA OH 1612 MHz survey 

OH/IR stars; v,<14 km/s (old pop.) 



(j){L)dL cx Lq ^■'^'^ ^ 



(18) 



that is, a luminosity function 4>{L) cx L^^-^'^ . This relation 
is fitted for Lq > O.375i?0(kpc)^ Jy kpc^, the detection 
limit. The exponent, a = —2.2, is close to the value of 
—2.0 obtained also in the Galactic center by Herman & 
Habing (1985). Chengalur et al. (1993) obtained a value of 
— 1.5 and Baud et al. (1981) obtained a value of —1.65, but 
their samples also contained young population (high tig), 
so the difference could stem from this. We are interested 
here only in the old population (older than 1.5 Gyr). A 
power law is generally accepted as a fit to the luminosity 
function, at least in the observed range. 

There is a cut-off for the faintest luminosities (Baud 
et al. 1981), which we take to be at Lq < O.2i?0(kpc)^ 
Jy kpc^, since the luminosity function peaks around 
Lq = O.2i?0(kpc)^ Jy kpc^ (Sjouwerman et al. 1998), and 
the number of fainter sources is lower (so the factor 
will make the number of nearby faint sources negligible) . 
Debattista et al. (2002) have also verified that the number 
of low luminosity sources is very small (corresponding to 
low luminosity is correlated with but the relation 
between both parameters cannot be used as a distance 
indicator since the dispersion of luminosities for a given 
Ve is large). Fig. 9 in Debattista ct al. (2002) shows that 
there are practically no low-luminosity stars (we < 6 km 
s^^; roughly less than O.2i?0(kpc)^ Jy kpc^ in average) at 
20° < 1^1 < 45° in the ATCA/VLA survey. Also, Sevenster 
et al. (1997a, fig. 6b; 2001, fig. 5b) find a very low number 
of low-Ue sources in the whole survey. 

With the luminosity function and the densities of the 
L02 or B02 models, we can calculate the number of counts 
in a given region of the sky as a function of R, and com- 
pare them with the distribution obtained from the data 
through cq. H16I) . Fig. E] shows the result for i? > 4 kpc 
in 12.8° < 1^1 < 45°. This figure shows an agreement be- 
tween models and data at i? > 4 kpc which corroborates 
that the luminosity function is approximately correct. 
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Fig. 14. Luminosity distribution of ATCA/VLA-1612 
MHz survey sources in the centre of the Galaxy < 
3°,|6| < 3°) and with radial velocities \vr\ > 30 km s^^. 
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Fig. 15. Distribution of ATCA/VLA-1612 MHz survey 
sources (fiux larger than 0.375 Jy, outfiow velocities < 
14 km s-i) in the region 12.8° < |;| < 45°, 0.75° < \b\ < 
2.5° as a function of the galactocentric radii for i? > 4 kpc. 
Model normalized to give the same number of sources than 
the data. Error bars stand only for Poisson errors. There 
is a good agreement with the B02 or L02 models (it does 
not matter whether it has a flat distribution or not, since 
this does not affect the regions at ii > 4 kpc) in this range 
of R. 



5.3. Deficit of stars in the inner disc 

If we apply this technique to the inner disc regions, we 
can see whether the model of ||31does or does not fit bet- 
ter than the L02 or B02 exponential models. We select the 
region \l\ > 12.8° (i.e. R > 1.75 kpc), 0.75° < |6| < 2.5° 
because this is far enough away from the bulge; it crosses 
the inner disc and is far enough away from the plane 
to consider the contribution of other possible plane old- 
population components (bar, ring) small in comparison 
to the error bars of the data. Rather than a plot as a 
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Fig. 16. Distribution of ATCA/VLA-1612 MHz survey 
sources (flux greater than 0.375 Jy, outflow velocities 
Ve < 14 km s~^) in the region 12.8° < |^| < /max, 
0.75° < l&l < 2.5°. Note that the points are not inde- 
pendent, however the point at Imax = 18° is enough to 
reject the B02 and L02 models beyond the 5a level. 



function of R (as in Fig. I15|) . we will show the ratio be- 
tween lower and higher R sources, between the sources 
in the inner disc {R < 3.75 kpc) and the remainder. 
The Poisson errors are very large if we divide the data 
into bins with different R. By means of eq. ()16|l . the ra- 

,. ( jV(fl<3.75 kpc) \ _ / jV([t,./sini]>190 km/s) \ . , . 
\N{B.>3.75 kpc) ) - {N{[vjsinl]<190 km/s) j SflOWn m 

Tablc[21 The lower the maximum longitude, the more sen- 
sitive we are to the inner disc population; however, the 
Poisson errors are larger (because we have fewer sources). 
Restricted to |/| < 18°, we get the larger discrepancy: 
the differences of an exponential model and the results 
from the data are around 5cr, significant enough. Figure 
^1 shows this graphically. 

The model of this paper gives a much better fit: all 
the ranges in Table [3 are within the 1.8a level. Perhaps 
the deficit of stars should be larger than that calculated 
in ^ but since we have not taken into account the errors 
of eq. H13() and the interpolation of the transition region, 
the compatibility with the flat distribution described in 
this paper is not excluded. 

Therefore, the conclusion of this section is again that 
a deficit of stars is necessary in the old population (older 
than 1.5 Gyr) of the disc with a density amplitude similar 
to that derived in fJ2| or even lower. It is noteworthy that 
this deficit of old population OH/IR stars has already been 
argued by Baud et al. (1981) for the populations around 1 
Gyr old. We quote a paragraph from §6.4 of this papcr^: 

That the sources with large AV exhibit the same 
kind of overall density distribution is to be ex- 
pected when they are as young as we have esti- 
mated. But the similarity of the density distribu- 



^ Baud et al. call the separation between the two peaks of 
OH/IR sources at 1612 MHz AV = 2v^ and use Rq = 10 kpc. 



tion of the sources with small AV(w 10^ yr old) 
and the extreme Population I objects is surpris- 
ing. We interpret this to mean that the large-scale 
CO distribution in the disc and hence the region 
of active star formation has not changed signifi- 
cantly during the last 10^ yr, or 4-8 rotations of 
the Galaxy. Therefore the 'hole' in the gas and in 
the OH/IR sources distribution at i? < 4.5 kpc and 
the maximum density at i? = 5 kpc appear to be 
rather stable phenomena. 

We corroborate these words for even older stars. 



6. Counts with bright sources in the mid-infrared 

In this section, we compare the disc models with bright 
mid-infrared star counts in the plane (i.e. a young popu- 
lation). In the mid- infrared, there is practically zero ex- 
tinction (^70 smaller than in V) even in plane regions. 
This means that extinction in the plane (b=0°) is typi- 
cally 0.1-0.2 mag in the mid-infrared, except in the very 
centre, where it can reach 0.4 mag. 

The comparison is carried out in Fig. 1171 In this case, 
MSX counts at 14.6 fim up to magnitude 3.0 are com- 
pared with different models. This is compared with the 
B02 model prediction (the luminosity function and the 
extinction were taken from the SKY model; updated ver- 
sion of Wainscoat et al. 1992) and with the extrapolation 
of the outer disc according to L02. Obviously, the L02 
and B02 models fail to reproduce the data for |Z| < 30°. 
COBE/DIRBE data at 12 /zm, representative of the flux 
of the young population, also give a nearly flat light distri- 
bution at 10° < |/| < 30°, 6 = 0° (Hammersley et al. 1994, 
fig. 1) instead of increasing towards / = 0°. Comparison 
with the models shows that the inner truncation of the 
exponential increase in the density in the disc can explain 
the difference. The model of this paper, which follows eq. 
(|13|l for i? < 4 kpc, can roughly fit the counts. Errors in 
the extinction are very low and cannot explain the ob- 
served difference between data and exponential models. 
The bulge, spiral arms, ring and in-plane bar also have 
some contribution which is not calculated here; note, how- 
ever, that they would increase the number of counts but 
never decrease it, so the interpretation of a deficit of stars 
is unavoidable. No attempt is made to fit the counts for 
other latitudes because we would need the scaleheight and 
the flare of the observed population with MSX-14.65 /im, 
dominated by a young population, and this is at present 
unknown. 



7. Summary of results: a smooth model for outer 
and inner disc 

We have used three different methods: 1) the inversion of 
the red clump giant distribution in near-infrared colour- 
magnitude diagrams to obtain the star density along the 
line of sight; 2) the determination of the density distri- 
bution of 1612 MHz sources by means of the distance 
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Table 2. The [^j^^^-^j ratio in the ATCA/VLA 1612 MHz survey and the predictions of three different models. 

Range of latitudes: 0.75° < \b\ < 2.5°. Errors in the data include Poisson errors and the errors due to the dispersion 
of values in R calculated according to eq. (|17|l [estimated by means of a Monte Carlo simulation] . 



Long, range Data Model "This paper" Model L02 Model B02 



12.8° < \l\ < 45° 0.65±0.22 0.46 0.81 0.80 

12.8° < \l\ < 30° 1.00±0.43 0.99 1.73 1.69 

12.8° < \l\ < 25° 1.00±0.48 1.59 2.91 2.79 

12.8° < |;| < 20° 1.10±0.62 1.97 3.83 3.71 

12.8° < |;| < 18° 1.00±0.64 2.10 4.30 4.11 

12.8° < |;| < 15° 1.00±0.81 2.28 4.84 4.70 



-0.5°<b<0.5° 

300 I — ' — ' — , — ' — ' — , — ' — / I V 




-45 -30 -15 15 30 45 

I (deg.) 

Fig. 17. MSX star counts in the plane at 14.6 fj,m up to 
magnitude 3.0. Comparison with the B02 disc model and 
the extrapolation of the L02 model of outer disc towards 
the centre with an exponential and flat density distribu- 
tion (this paper). All models were normalized with the 
counts to fit the counts at I — 45°. It is clear that the 
deficit of stars in the young stellar population represented 
here is necessary. Note that the models represent only the 
disc (without a bulge, bar or ring) while the counts include 
everything in the line of sight. 




Fig. 18. Contour diagram of logj^Q p (pc^^) in the yz-plane 
of the Galaxy (perpendicular to the line Sun-Galactic cen- 
tre) with y between —12.0 and 12.0 kpc and z between 
— 1.0 and 1.0 kpc (vertical scale in the plot multiphed by 
a factor 5). Lower contour: log^g P (pc~^)=-2.1; step=0.15. 
Model: L02 for i? > 4 kpc, and eq. JTSJ for i? < 4 kpc. 
The abrupt transition (at \y\ = R = 4 kpc) is due to the 
discontinuous change of regime between the inner disc and 
the outer disc. Note also that the model of this paper is 
valid only for regions not very far from the plane. 



determination of OH/IR sources from their kinematical 
information; 3) the analysis of near-infrared star counts 
and comparison with models. All the tests agree with the 
following picture for the inner disc: 

— There is a deficit of stars with respect to an exponential 
disc in the Milky Way's inner stellar disc for 2.25 kpc 
< i? < 4 kpc. The density in the plane {b « 0°, i.e. 
z ~ 0) is almost independent of R instead of being an 
exponential law of the type p oc exp(— i?//i). 

— There is a slight increase in scaleheight in the central 
disc, hz = 509 — 48i?(kpc) pc, so the combined 
effect of deficit of stars in the plane and higher 
scaleheight produces an important deficit of stars for 
low latitude regions, although not for high latitude 
regions. Roughly, the density follows an expression of 
the type p{R,z) w 5.1 ± O.6jO0e™9±'5^-*s±2OH(fcpc) pc ^ 



a result that is also compatible with a flared inner 
disc with a fixed scalelength of 2.4 kpc: p{R, z) « 

/)„p 2.4kpc p 509±67-48±20-R(fcpc) pc P'" 

^ 509±67-48±20fl(fcpc) pe- 

lt is our objective to investigate further the real 
nature of the inner stellar disc, with a more complete 
analysis of 2MASS All Sky Release star counts. See the 
graphical representation in Fig. ^| (both possibilities 
are numerically equivalent). Note that these numbers 
are also subject to several systematic errors (^ 15% in 
the density between latitudes 2 and 6 degrees), so the 
parameters might change; although not as much as to 
allow the fit of a purely exponential disc (which would 
require errors greater than 70%). Moreover, the deficit 
of stars is corroborated with the OH/IR sources and 
the star counts. These parameters are applicable only 
for regions not very far from the plane (|6| < 6°). 
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Fig. 19. Contour diagram of log^Q p (pc^'^) of a possible 
interpolated/extrapolated smooth model of the disc ac- 
cording to eq. (Uni) for R < Rq and L02 for R > Rq 
in the j/z-plane of the Galaxy (perpendicular to the line 
Sun-Galactic Centre) with y between —12.0 and 12.0 kpc 
and z between —1.0 and 1.0 kpc (vertical scale in the plot 
multiplied by a factor 5). Lower contour: log^o p (pc~"^) = 
-2.1; step = 0.15. 

— Both the young and old population have this deficit 
of stars, and at present no significant differences in 
their distributions have been deduced from the data. 
Presumably, a more accurate measure of the distribu- 
tion might reveal some differences. At present we can 
say that, roughly speaking, the deficit of stars with re- 
spect the extrapolation of an exponential distribution 
is probably a rather stable feature of the disc, which 
might be due to the existence of an in-plane bar that 
sweeps out the near-plane stars. 

— An expression that summarizes the L02 outer disc 
(i? > 6 kpc) and the inner disc (2.5 kpc < i? < 4 
kpc) with a smooth transition between two regimes is: 



p{R,z) 



PqC' 



X e 



(19) 



with 

h,XR) « 285[1 + 0.21 kpc~^(i?-i?0) 
4- 0.056 kpc"^(i? - Rq)^] pc. 



(20) 



This density is plotted in Fig.^| Although this expres- 
sion was not tested for all ranges (for instance, it is not 
tested around R — 5 kpc or for i? < 2 kpc), this can 
be used as an approximate interpolation/extrapolation 
that models the disc in i? < Rq. For the outermost 
disc, it is better still to use the expression of L02. The 
innermost disc {R < 2 kpc) reveals with the expres- 
sion (|19|l a very deep hole. Unfortunately, we cannot 
be sure of the existence of this deep hole there, since 
the disc and the bulge cannot be easily separated in 



that region. It is clear at least that between 2 and 4 
kpc the deficit of stars is very significant. 
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